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INTRODUCTION 


The  influences  of  curvature  on  flow  behavior  are  important  to  a 
variety  of  machines  including  pump  impellers,  radial  compressors,  nozzles 
and  diffusers,  axial  flow  compressors,  and  axial  flow  turbines.  To 
develop  predictions  of  these  flows,  bench  mark  experimental  results  are 
required.  The  present  study  is  part  of  an  experimental  effort  to  achieve 
this  aim.  In  it,  flow  instabilities  in  a  constant  cross  section  channel 
with  constant  curvature  are  examined.  These  results  not  only  enhance 
understanding  of  flows  in  more  complex  geometries,  but  also  provide  a 
building  block  for  the  development  of  improved  prediction  schemes  and  more 
efficient  machine  designs. 

A.  BACKGROUND 

Until  recently,  only  two  laminar  flows  were  known  to  occur  prior  to 
the  onset  of  turbulence  in  curved  channels.  At  low  enough  Re,  the  most 
basic  laminar  flow  occurs.  This  is  curved  channel  Poiseuille  flow  (CCPF). 
At  higher  Re,  secondary  flow  develops  due  to  centrifugal  instability  of 
CCPF.  This  second  laminar  flow  consists  of  an  array  of  streamwise- 
oriented,  counter-rotating  pairs  of  vortices  and  is  called  Dean  vortex 
flow. 

If  the  spanwise  spacing  of  vortex  pairs  is  X,  then  a  non-dimensional 
spanwise  wavenumber  if  the  vortex  pairs  is  a=2  (d/2)7r/X.  For  given  a  and 
radius  ratio  q  ,  classical  linear  stability  analysis  predicts  the  neutrally 
stable  De  which  sets  the  lower  limit  where  such  vortices  may  be  present. 
The  minimum  of  the  neutral  stability  curve  for  the  primary  instability 
then  occurs  at  Dec  and  ac.  For  q =0.979,  Dee=36,92. 

Recent  numerical  work  by  Finlay  et  al.  [Ref.  1]  predicts  these  Dean 
vortices  become  unstable  at  higher  Re  to  two  types  of  travelling  waves, 
resulting  in  two  kinds  of  wavy  vortex  flows,  referred  to  as  undulating 
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Dean  vortex  flow  and  twisting  Dean  vortex  flow.  The  first  mention  of 
experimentally  observed  waviness  is  given  by  Kelleher  et  al.  [Ref.  2],  and 
the  first  photographic  evidence  of  vortex  motions  believed  to  be  twisting 
is  presented  by  Ligrani  and  Niver  [Ref.  3].  Additional  experimental 
evidence  of  wavy  vortex  flow  in  a  rotating,  curved  channel  is  given  by 
Matsson  and  Alfredsson  [Ref.  4].  Finlay  and  Nandakumar  [Ref.  5]  provide 
experimental  evidence  of  undulating  Dean  vortex  flow  and  twisting  Dean 
vortex  flow  which  are  compared  to  wavy  vortex  results  obtained  by  direct 
simulation  of  the  time-independent  three-dimensional  Navier-Stokes 
equations.  The  experimental  results  were  obtained  using  the  same  large 
aspect  ratio  curved  channel  described  by  Ligrani  and  Niver  [Ref.  3], 

Another  category  of  motions  which  has  been  observed  is  related  to 
appearances  and  disappearances  of  Dean  vortex  pairs  across  the  span  of  the 
channel.  Evidence  of  such  behavior  for  channel  flow  was  first  provided  by 
Ligrani  and  Niver  [Ref.  3]  in  the  form  of  photographic  sequences  showing 
small  secondary  vortex  pairs  observed  temporarily  in  spanwise/radial 
planes  for  De  ranging  from  75  to  220.  Additional  study  of  different 
phenomena  that  result  in  these  appearances  and  disappearances  is  presented 
by  Ligrani  et  al.  [Ref.  6J. 

B.  OBJECTIVES 

The  present  study  focuses  on  the  behavior  of  Dean  vortices  .in 
rectangular  curved  channels  with  40  to  1  aspect  ratios.  Methods  are 
employed  to  study  flows  at  De  from  49.2  to  450.  This  range  of  De 
represents  flow  regimes  from  laminar  CCPF  to  fully  turbulent  flow.  Of 
particular  interest  are  undulation,  twisting,  appearances  of  vortex  pairs, 
disappearances  of  vortex  pairs,  and  their  connections  to  the  development 
of  fully  turbulent  flow. 

C.  ORGANIZATION 

Chapter  II  describes  the  construction  and  operation  of  the  curbed 
channels  used  for  the  present  study.  Included  is  information  on  the 
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determination  of  bulk  flow  rates  in  the  channels.  Chapter  III  describes 
equipment  used  for  measurement  of  pressure  and  velocity  surveys  plus 
techniques  used  to  obtain  real  time  flow  visualization  results.  Chapter 
IV  presents  experimental  results.  Chapter  V  presents  the  summary  and 
conclusions.  Appendix  A  contains  uncertainty  analysis  results.  Appendix 
B  describes  the  software  employed.  Appendix  C  contains  sequential  Dynamic 
Signal  Analyzer  steps  to  obtain  power  spectra  plots.  Appendix  D  contains 
all  tables.  Appendix  E  contains  photographs,  schematics  and  result 
figures.  Appendix  F  contains  individual  power  spectra  plots. 
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IZ.  EXPERIMENTAL  FACILITIES 


A.  CURVED  CHANNEL  DESCRIPTIONS 

Two  nearly  identical  curved  channels,  a  transparent  channel  and  a  heat 
transfer  channel,  were  utilized.  A  photograph  of  the  transparent  channel, 
a  photograph  of  the  blower  assembly,  and  a  schematic  of  the  channel  are 
shown  in  Figures  1,  2  and  3,  respectively.  The  channel  is  described  in 
detail  by  Ligrani  and  Niver  [Ref.  3],  The  facility  is  an  open-circuit 
suction  facility  designed  for  low-speed  transition  studies.  At  the  inlet, 
an  aluminum  honeycomb  and  three  screens  reduce  spatial  non-uniformities  in 
the  flow.  These  are  followed  by  a  two-dimensional  nozzle  with  a  20  to  1 
contraction  ratio  which  accelerates  the  flow  into  a  2.44  m  long  straight 
duct.  The  straight  section  provides  fully  developed  laminar  flow  at  the 
inlet  of  the  curved  section  which  then  follows,  for  Dean  numbers  up  to  640 
[Ref.  3).  The  radii  of  the  concentric  convex  and  concave  walls  in  the 
curved  section  are  596.9  mm  and  609.6  mm,  respectively.  The  180  degree 
curved  section  is  then  followed  by  another  2.44  m  long  straight  duct,  and 
then  a  honeycomb,  several  screens,  a  diffuser  and  an  outlet  plenum.  From 
the  inlet  f'f  the  first  straight  section  to  the  outlet  of  the  second 
straight  section,  the  interior  dimensions  of  the  facility  are  1.27  cm  for 
the  height  and  50.80  cm  for  the  width,  providing  an  aspect  ratio  of  40  to 
1.  Tolerances  for  the  height  and  width  are  0.015  cm  and  0.050  cm, 
respectively,  maintained  using  longitudinal  and  cross-beam  supports  along 
the  length  of  the  test  section.  These  supports  are  com.  cted  to  two 
continuous,  seamless  sheets  of  polycarbonate  (commercial  name  is  Lexan) 
which  form  the  two  primary  interior  walls  adjacent  to  the  flow  stream. 
The  outlet  plenum  (labeled  #1  in  Figure  3)  is  connected  to  a  second  plenum 
(labeled  #2  in  Figure  3)  using  50.8  mm  diameter  piping.  Plenum  #2  is 
maintained  at  low  pressure  by  a  1/3  horsepower  (250  W)  blower  manufactured 
by  ICG  Industries,  Inc.  A  brass  gate  valve  is  used  to  adjust  the  flow 
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rate,  which  is  determined  from  the  pressure  drop  across  an  orifice  plate 
located  in  the  pipe  between  the  two  plenum  chambers. 

In  the  transparent  channel  a  slot,  3.2  mm  wide  and  76.2  mm  long  is 
located  in  the  convex  wall.  It  is  aligned  in  the  spanwise  direction  50.8 
mm  off  the  centerline  at  a  location  120  degrees  from  the  start  of 
curvature.  A  support  block  of  polycarbonate  (152.4  mm  by  63.5  mm  by  22.2 
mm)  is  joined  using  epoxy  to  the  wall  surrounding  the  slot  to  maintain  the 
convex  wall's  dimensional  integrity,  and  to  minimize  the  deflection  of  the 
curved  wall  in  the  vicinity  of  the  slot.  Foam  is  used  to  line  the  slot  to 
allow  probe  insertion  with  no  air  leakage  to  the  channel  interior.  This 
foam  is  carefully  positioned  to  avoid  any  bumps  or  discontinuities  in  the 
flow.  Additional  details  are  given  by  Ligrani,  et  al.  (Ref.  7] 

In  the  second  curved  channel,  a  slot  identical  to  the  one  previously 
described  is  located  straddling  the  centerline,  aligned  in  the  spanwise 
direction,  and  positioned  19.05  cm  downstream  of  the  convex  surface.  This 
channel  is  configured  for  heat  transfer  measurements,  but  the  Blot 
location  is  the  only  significant  difference  pertaining  to  this  study. 
Additional  details  are  given  by  Hughes.  [Ref.  8) 

B.  DEAN  NUMBER  DETERMINATION 

Dean  number,  channel  Reynolds  number,  channel  bulk  velocity  and  mass 
flow  rate  are  determined  by  measuring  the  pressure  drop  across  the  3.81  cm 
(1,5  inch)  diameter  ASME  sharp  edged  orifice  plate.  Also  measured  are 
barometric  pressure,  and  ambient  temperature.  The  expansion  coefficient 
Y  is  then  calculated  from  Holman  and  Gajda  [Ref.  10).  An  initial  Dean 
number  guess  is  then  made  by  interpolating  a  table  of  Dean  number  versus 
orifice  pressure  orop  compiled  by  Niver  [Ref.  11).  The  density  of  air  is 
computed  using  the  ideal  gas  law.  An  initial  mass  flow  rate  is 
subsequently  calculated  based  on  the  initial  Dean  number.  The  Reynolds 
number  through  the  pipe  is  calculated  using  the  equation  given  by: 


5 


i?e  =/ft— %- 
p  pv^E 


Flow  coefficient  K  is  interpolated  from  ASME  tables  [Ref.  12].  With  this, 
mass  flow  rate  is  determined  using  an  equation  of  the  form: 


lh=KA0Yj2pKP7r 

Iterations  continue  until  convergence  to  within  1%  of  previous  mass  flow 
rate.  Channel  bulk  velocity,  channel  Reynolds  number  and  Dean  number  are 
then  calculated  using  the  equations: 


U= 


rh 

PAch 


and  De=Rech, 

respectively.  A  plot  of  pressure  drop  versus  De  is  shown  in  Figure  4  for 
a  barometric  pressure  of  29.92  inches  and  an  ambient  temperature  of  19 
degrees  Celsius. 
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III.  EXPERIMENTAL  PROCEDURES 


A.  STREAMWISE,  SPANWISE  AND  RADIAL  VELOCITY  SURVEYS 
1.  Five-hole  Pressure  Probe  Description 

A  photograph  and  two  schematic  diagrams  of  the  miniature  five-hole 
pressure  probe  are  shown  in  Figures  5,  6  and  7,  respectively.  The 
diameter  of  the  probe  tip  io  1.38  mm.  The  probe  is  arranged  such  that  the 
central  tube  is  surrounded  by  four  other  tubes.  The  ends  of  each  of  the 
surrounding  tubes  are  beveled  at  a  45°  angle  with  respect  to  the  mouth  of 
the  central  tube.  Each  of  the  five-hole  tubes  is  made  from  stainless 
steel  hypodermic  tubing  with  inner  and  outer  diameters  of  0.203  mm  and 
0.406  mm,  respectively.  The  probe  was  constructed  by  first  spot  welding 
together  these  five  tubes.  The  tubes  were  then  bent  in  a  90°  turn  with  a 
radius  of  1.52  mm  to  avoid  any  tube  crimping.  Silver  solder  was  then  used 
to  join  the  five  tubes  together  continuously  along  their  length.  Finally, 
portions  of  the  probe  tip  were  beveled  to  appropriate  angles  using  a 
surface  grinder.  This  approach  is  slightly  different  than  the  one 
described  by  Treaster  and  Houtz.  [Ref.  13] 

Figure  6  shows  that  the  present  probe  tip  extends  5.08  mm  (4  tip 
diameters)  beyond  the  bend  in  the  stem  to  minimize  the  effects  of  stem 
flow  blockage  at  the  measuring  location.  Similarly,  the  tip  is  22.9  mm 
(19  tip  diameters)  from  the  9.53  mm  diameter  mounting  cylinder.  Within 
this  cylinder,  each  tube  is  connected  to  a  1.60  mm  diameter  tube  using  a 
small  manifold  for  each  tube  connection.  The  mounting  cylinder  was  sized 
for  easy  handling  and  mounting.  On  top  of  the  cylinder  is  a  second  set  of 
manifolds  with  five  male  connectors  for  attachment  to  6.35  mm  outer 
diameter/4.67  mm  inner  diameter  Imperial  Eastman  Polyflo  tubing. 

The  miniature  five-hole  probe  was  designed  to  minimize  the  adverse 
effects  of  flow  blockage  and  spatial  resolution  when  used  for  measurements 
in  the  curved  channel.  When  in  place  for  measurements  at  the  centerline 
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of  the  channel,  the  mounting  cylinder  is  not  exposed  to  the  flow,  and  the 
probe  produces  4.8%  blockage  over  a  12.7  mm  x  12.7  mm  area.  A  five-hole 
rather  than  a  four-hole  design  was  chosen  since  it  was  easier  to  calibrate 
and  because  corrections  for  spatial  resolution  are  more  readily  applied 
[Ref.  14].  The  tip  was  conical  rather  than  prismatic  to  minimize  near¬ 
wall  flow  blockage  effects,  which  are  particularly  severe  for  prismatic 
type  probes. 

2.  Experimental  Details 

Figure  8  shows  a  schematic  of  the  measurement  chain  used  for  the 
five-hole  pressure  probe.  Each  tube  of  the  five-hole  probe  iB  connected 
via  flexible  plastic  tubing  (less  than  20  cm  in  length  to  minimize  the 
volume  of  air  in  the  tubes)  to  a  Validyne  model  DP103-06  variable 
reluctance  differential  pressure  transducer.  A  separate  transducer  is 
used  to  measure  pressure  relative  to  the  atmosphere  for  each  port  of  the 
five-hole  probe.  Each  transducer  has  a  full-scale  range  of  2.5  mm  of 
water  differential  pressure.  Shielded  cables  connect  the  transducers  to 
separate  Celesco  CD01D  carrier  demodulator  units.  The  carrier  demodulator 
units  produce  DC  voltage  output  proportional  to  the  pressure  signals  from 
the  transducers.  Each  carrier  demodulator  output  is  connected  to  a 
Hewlett  Packard  (HP)  3498A  extender  which  is  controlled  by  an  HP  3497A 
data  acquisition/control  unit.  The  HP  3497A  is  controlled  by  an  HP  Series 
300,  Model  236  computer  which  also  processes  the  data.  Voltage  samples 
are  recorded  sequentially  from  the  five  channels  using  an  overall  sampling 
rate  of  about  10  Hz. 

3 .  Probe  Calibration 

The  probe  was  calibrated  by  Fields  [Ref. 15]  using  a  procedure 
similar  to  Ligrani,  et  al.  [Ref.  7].  During  calibration,  the  probe  was 
mounted  in  a  calibration  sled  and  placed  in  an  open-circuit  blower  tunnel 
which  provided  uniform  flow  in  the  test  section.  A  manually  operated 
traverse  was  used  to  position  the  probe  through  a  five  yaw  angle  and  five 
pitch  angle  array.  Total  pressure  was  measured  with  a  Kiel  probe,  and 
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tunnel  static  pressure  was  measured  near  the  probe  position  using  wall 
static  taps.  For  each  probe  position  during  calibration,  yaw,  pitch, 
total  and  total  minus  static  pressure  coefficients  were  calculated. 

4.  Curved  Channel  Procedures 

a.  Traversing  Mechanism 

When  measurements  are  made  in  the  channel,  the  probe  is 
oriented  using  a  mount  which  allows  initial  yaw  and  pitch  angles  to  be 
adjusted  within  ±0.1°  and  ±0.5°,  respectively.  The  mounting  block  is 
attached  to  a  travelling  block  which  is  part  of  an  automated  two- 
dimensional  (X-Y)  traversing  mechanism.  Moving  parts  are  driven  by  20 
thread  per  inch  (0.8  thread  per  mm)  drive  screws  which  are  rotated  by 
Superior  Electric  type  M092-FD310  stepping  motors.  These  motors  are 
controlled  by  a  MITAS  Motion  Controller,  which  directs  movement  of  the 
probe  according  to  instructions  from  the  computer. 

b.  Measurements 

The  probe  is  positioned  and  aligned  before  measurements  are 
made.  The  probe  is  first  inserted  in  the  slot  of  the  insulated  channel  at 
a  spanwise  location  25.4  mm  from  the  centerline  so  that  it  just  touches 
the  straight  wall  farthest  from  the  slot.  It  is  then  visually  aligned 
with  the  flow  and  then  moved  to  the  channel  position  at  the  midpoint  of 
the  spanwise-radial  plane  using  the  MITAS  controller.  Small  adjustments 
are  made  to  the  yaw  and  pitch  angle  by  adjusting  the  probe  mount. 
Software  program  ORIENT,  described  in  Appendix  B,  is  used  to  check  yaw  and 
pitch  coefficients  as  this  is  undertaken  until  the  yaw  coefficient  is  zero 
and  the  pitch  coefficient  is  sufficiently  small.  The  final  orientation  of 
the  probe  gives  a  yaw  coefficient  of  0.0±0.01  and  a  pitch  coefficient  of  - 
0.4.  This  is  equivalent  to  +3.0°  offset.  More  accurate  pitch  alignment 
forces  the  probe  against  the  edge  of  the  slot,  hence  the  pitch  offset  is 
accounted  for  during  data  reduction.  No  offset  is  applied  or  needed  for 
the  yaw  angle. 


9 


In  a  typical  spanwise/radial  traverse,  the  probe  moves  in  the 
radial  direction  from  the  concave  wall  to  the  convex  wall  as  measurements 
are  made.  To  begin  measurements  at  a  new  spanwise  location,  the  probe 
first  moves  2.54  mm  in  the  spanwise  direction,  then  in  the  reverse  span 
direction  1.27  mm,  and  finally  toward  the  concave  wall  until  it  is  1.27  mm 
away.  A  new  radial  or  spanwise  measuring  traverse  then  begins.  This 
procedure  reduces  the  adverse  effects  of  backlash  and  gives  repeatable 
data  regardless  of  the  direction  of  spanwise  probe  motion. 

When  measurements  are  made  at  a  particular  location,  50 
samples  of  each  pressure  are  averaged.  This  is  done  by  sampling  pressure 
from  probe  ports  in  sequence  so  that  a  long  averaging  time  is  provided  for 
each  channel,  thereby  minimizing  the  influence  of  brief  unsteadiness  in 
the  channel.  For  each  probe  location,  about  3  min  of  acquisition  and 
processing  time  is  required.  A  traverse  over  an  entire  plane  requires 
about  8  h.  The  entire  process  is  computer  controlled  and  generally 
performed  at  night  when  there  is  minimal  unsteadiness  at  the  channel  inlet 
from  laboratory  disturbances, 
c.  Data  Reduction 

Measured  pressures  from  each  probe  port  are  stored  by  the 
computer  during  each  traverse.  These  data  are  then  corrected  to  account 
for  the  effect  of  finite  spatial  resolution  using  procedures  described  by 
Ligrani  et  al.  [Ref.  14].  The  procedure  corrects  pressures  to  account  for 
the  fact  that  they  are  not  measured  at  the  same  physical  location  by 
adjusting  pressures  from  ports  2,  3,  4  and  5  so  that  they  appear  to  be 
measured  at  port  1.  The  correction  primarily  affects  secondary  flow 
velocity  vectors  in  the  plane  normal  to  the  bulk  flow  direction.  With  the 
corrections,  these  secondary  vector  magnitudes  may  change  as  much  as  15% 
at  locations  where  velocity  gradients  are  largest. 

Yaw  and  pitch  pressure  coefficients  are  then  calculated  for 
each  location.  Using  these  coefficients  and  measured  pressures,  pitch  and 
yaw  angles  are  calculated  at  each  location.  The  total  minus  static 
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pressure  coefficient  is  then  calculated.  From  this  total  velocity 
magnitude  is  calculated  using  the  equation: 

U= 


The  three  velocity  components  are  subsequently  given  by: 

U0=Ucos(y)cos(p) 

Ur=Usin(y) 

and  Uz=Ucos(y)sin(p) 

where  y  and  p  are  the  yaw  and  pitch  angles,  respectively.  Mean  vorticity 
components  are  then  calculated  using  finite  difference  schemes  applied  to 
the  three  components  of  mean  velocity.  Equations  for  the  vorticity 
components  as  functions  of  velocity  are  given  by: 


we 


(i) 


z 


and 


2  i 


Hi 

dr 


Streamwise  derivatives  are  neglected  when  calculating  radial  and  spanwise 
vorticities.  In  the  final  phase  of  data  reduction,  corrections  are  made 
to  account  for  deflection  of  flow  streamlines  near  the  probe  tip  in  the 
presence  of  transverse  gradients  of  the  streamwise  velocity.  This  is 
referred  to  as  the  downwash  velocity  correction  [Ref.  14],  and  is  based  on 
the  idea  that  the  downwash  velocity  in  a  particular  direction  scales  with 
-he  transverse  gradient  of  streamwise  velocity  in  the  same  direction. 
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B.  CONVEX  NALL  TURBULENCE  SURVEYS 

1.  12  Element  Hot-film  Probe  Description 

A  photograph  of  the  hot-film  array  is  shown  in  Figure  9.  The 
array  has  12  hot-film  anemometers  mounted  equidistant  from  each  other  on 
a  specially  designed  aluminum  frame.  Probe  locations  are  tabulated  in 
Table  1,  which  indicates  that  probe  1  is  positioned  56.95  mm  from  the 
channel  spanwise  centerline.  Each  probe  is  silver  plated  platinum  with  a 
diameter  of  0.625  microns  and  an  active  length  of  0.102  mm.  The  frame  is 
designed  to  fit  snugly  into  the  slot  located  120°  from  the  start  of 
curvature  in  the  transparent  curved  channel.  With  thin  shims  placed  under 
the  ends  of  the  frame,  the  probe  is  mounted  so  that  the  measuring  plane  of 
the  gages  is  flush  with  the  convex  wall. 

2.  Experimental  Details 

Figure  10  shows  a  schematic  of  the  measurement  chain,  including 
the  probes,  and  data  acquisition  and  processing  equipment.  With  the 
present  setup,  any  six  probes  can  be  connected  to  each  of  six  different 
DANTEC  56C17  Constant  Temperature  Anemometer  bridge  units  shown  in  the 
photograph  Figure  11.  Bridge  outputs  are  sent  simultaneously  to  a  HP 
3497A  Data  Acquisition/Control  Unit  and  to  DANTEC  56N20  Signal 
Conditioning  Units,  in  the  discussion  which  follows,  the  former  signals 
are  referred  to  as  raw  signals,  whereas  the  latter  are  referred  to  as 
conditioned  signals.  The  output  of  each  signal  conditioning  unit  is 
directed  to  the  HP  3497A  Data  Acquisition/Control  Unit  which  is  controlled 
by  an  HP  Series  300,  Model  236  Computer  which  also  processes  the  data. 
With  the  control  unit  and  computer,  voltage  samples  are  recorded 
sequentially  from  the  six  channels  using  an  overall  sampling  rate  of  about 
5.9  Hz. 

3.  Curved  Channel  Procedures 

a.  Measurements 

The  data  acquisition  system  is  arranged  so  that  the  signal 
from  the  hot-film  gage  closest  to  the  spanwise  centerline  is  sampled 
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first.  Signals  from  adjacent  gages  across  the  span  of  the  array  are  then 
sampled  in  sequence.  The  sequence  is  repeated  until  1000  sampled  voltages 
for  each  gage  are  obtained.  Each  bridge  uses  an  overheat  ratio  setting  of 
1.6.  Each  signal  conditioning  unit  is  set  with  filters  to  pass  signals 
with  frequencies  from  0.1  Hz  to  1.0  kHz.  A  gain  setting  of  500  is  also 
used. 

b.  Data  Reduction 

After  sampled  voltages  are  stored  in  the  computer,  they  are 
reduced  by  a  factor  of  500  to  account  for  the  amplifier  gain.  Mean  values 
of  the  raw  and  conditioned  voltage  signals  are  then  determined.  In 
addition,  rms  voltage  for  each  conditioned  channel  is  also  determined. 
Normalized  rms  voltage  for  each  channel  is  computed  as  the  ratio  of  rms 
voltage  to  the  raw  mean  voltage. 

Covariances  are  determined  for  the  conditioned  voltage  signals 
using  the  equations 


n  _ 

COV=—'YJ  ( Erefi-Eref )  {Eobs^Eobs) 
n  i=l 

Here,  n  is  the  number  of  samples,  ref  is  the  reference  channel  and  obs  is 
the  observed  channel.  This  is  equivalent  to: 

e'  ( x )  e1  (x+xO 

Two  different  probe  configurations  were  utilized.  First,  covariances  for 
probes  2,  4,  6,  8,  10  and  12  with  respect  to  probe  2  located  at  z/d=4.73 
were  determined  for  De  from  50.25  to  399.8.  This  sequence  is  termed  'in 
to  out'.  Covariances  for  probes  12,  10,  8,  6,  4  and  2  with  respect  to 
probe  12  located  at  z/d=7,23  were  determined  for  De=50.25  to  399.8.  This 
sequence  is  termed  'out  to  in'.  Second,  'in  to  out'  covariances  for 
probes  1  through  6  with  respect  to  probe  1  located  at  z/d=4.48  were 
determined  for  De  from  49.2  to  400.3.  'Out  to  in'  covariances  for  probes 
6,  5,  4,  3,  2  and  1  with  respect  to  probe  6  located  at  z/d=5.73  were 
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determined  for  De=49.2  to  400.3.  After  determining  covariances,  spatial 
correlations  are  determined  using  the  equation: 


COR= 


_ COV _ 

( rmsEref )  ( imsEobs ) 


This  equation  is  equivalent  to: 

e'  (x)  e1  (x+x'i 
]/e'  (x)2\Je'  (x+x')2 


Length  scale,  representing  the  area  under  a  curve  of  correlation  versus 
z/d  is  determined  by  the  equation: 

n 

LS=Y  (COR)  (A Z) 

3  =  1 


or  alternatively: 


Lx= J  CORdz 


Here,  AZ  is  the  spanwise  distance  between  active  probes.  The  summation 
was  used  to  calculate  length  scales  and  is  an  approximation  to  the  exact 
integral  method. 

C.  FLOW  VISUALIZATION 

In  order  to  obtain  information  on  secondary  flow  motions  resulting 
from  the  vortices,  air  in  the  transparent  channel  was  partially 
contaminated  with  smoke  and  then  photographed.  The  device  used  to 
generate  the  smoke  is  described  by  Ligrani  and  Niver.  [Ref.  3]  Smoke 
patterns  were  illuminated  in  spanwise/radial  planes  using  spotlights 
directed  at  spanwise  slits  through  black  paper  used  to  line  the  convex 
surface  exterior  of  the  channel.  A  Sony  DXC-M3  video  camera  with  a 
Fujinon-TVZ  1:1.7/10-140  nun  DCL-914BY  zoom  lens  was  used  to  record  the 
illuminated,  time-varying  patterns  as  shown  in  Figure  12.  The  camera  was 
connected  to  a  Sony  VO-6800  portable  videocassette  recorder  which  collects 
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images  at  rates  as  high  as  60  frames  per  second.  Recorded  images  are  then 
played  back  using  a  Sony  VO-5800  videocassette  recorder  connected  to  a 
SONY  PVM-1910  Trinitron  color  monitor.  By  employing  the  scan/stop  action 
feature  of  this  recorder,  photographs  of  instantaneous  images  using  the 
setup  shown  in  Figure  13  were  obtained  and  then  arranged  in  sequences  to 
show  time-varying  phenomena.  New  time-varying  sequences  are  presented  in 
Chapter  IV.  This  chapter  additionally  gives  results  from  analysis  of 
video  tapes  previously  obtained  by  Ligrani  et  al.  [Ref.  6].  Data  were 
analyzed  at  De  from  60  to  150  at  angles  from  the  start  of  curvature  from 
85°  to  135°.  Of  primary  interest  were  frequencies  of  different  appearance 
and  disappearance  events  described  by  Ligrani  et  al.  [Ref.  6]. 
Tabulations  of  results  are  given  in  Tables  2  and  3.  In  making  these 
tabulations,  vortex  pairs  which  moved  into  or  out  of  the  field  of  view 
during  the  time  sequences  were  ignored.  All  frequencies  were  also 
converted  to  non-dimensional  form  using  the  equation: 

f  nfd 
nd~  U 

These  are  also  tabulated  in  Table  4. 

The  two  new  time-varying  sequences  presented  here  were  obtained  using 
the  setup  of  Figure  13  using  a  Canon  A-l  SLR  camera  with  a  35-70  mm,  f3.5- 
22  lens.  Settings  of  55  mm,  f5.6  and  speed  1/8  second  for  Kodak  Tri-X 
pan,  400  ASA  film  were  utilized  following  the  recommendations  of  Longest 
[Ref.  16] . 

D.  DEAN  NUMBER  SURVEYS  OF  POWER  SPECTRA  AND 

LONGITUDINAL  VELOCITY  FLUCTUATIONS 

I.  Experimental  Details 

A  DANTEC-  55P14  single  hot-wire  probe  is  used  to  obtain  surveys  of 
power  spectra  and  longitudinal  velocity  fluctuations.  The  probe  sensor  is 
platinum  plated  tungsten  wire  with  a  diameter  of  5  microns  and  a  sensor 
length  of  1.25  mm.  Because  support  prongs  are  bent,  they  are 
approximately  parallel  to  the  flow  direction  at  the  sensor  location,  which 
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r  places  the  sensor  normal  to  the  bulk  flow  direction.  During  measurements, 

'  the  probe  is  rigidly  mounted  in  an  X-Y  traverse  such  that  the  sensor 

I  itself  is  oriented  in  the  spanwise  direction  normal  to  the  flow  at  a 

streamwise  location  120°  from  the  start  of  curvature.  Surveys  were 

j 

obtained  at  each  of  five  locations.  These  locations  are  listed  in  Table 

l 

5  and  shown  in  Figure  16  relative  to  vortex  pair  structures.  The  probe  is 
operated  at  an  overheat  ratio  of  1.8  using  a  DISA  55M10  Constant 
Temperature  Anemometer  bridge.  The  output  from  the  bridge  is  routed 
through  a  DANTEC  56N20  Signal  Conditioning  Unit.  Settings  on  the  signal 
conditioner  were  varied  with  Dean  number  range  as  shown  in  Table  6.  The 
output  from  the  signal  conditioner  was  routed  to  a  DANTEC  56N22  Mean  Value 
Unit  set  at  a  10  second  integration  interval.  The  output  from  the  signal 
conditioner  was  additionally  sent  to  a  Hewlett  Packard  6944A 
Multiprogrammer.  This  high  speed  data  acquisition  system  was  controlled 
by  a  Howlett  Packard  300  Series  Model  9000  computer  which  also  stored  and 
processed  the  data.  Acquisition  sampling  rate  was  2500  Hz  if  De  was  less 
than  275  or  5000  Hz  if  Dean  number  was  higher  than  275. 

Another  output  from  the  bridge  was  directed  to  a  second  signal 
conditioner  operated  with  amplifier  gain  of  500,  low  pass  filter  of  1.0 
kHz,  and  a  high  pass  filter  of  0.1  Hz.  The  output  from  this  signal 
conditioner  was  then  sent  to  a  DISA  56N25  RMS  Unit  set  with  a  10  second 
integration  interval. 

■  After  the  DISA  56N20  Signal  conditioner,  signals  were  also  sent  to 

J  the  Hewlett  Packard  3562A  Dynamic  Signal  Analyzer.  This  device  is  used  to 

i 

j  obtain  spectra  of  the  time-varying  signals  from  the  hot-wire  bridge. 

j  Figure  14  shows  the  setup  for  measurement,  data  acquisition  and 

processing.  Figure  15  gives  an  overview  of  the  internal  measurement  chain 
of  the  dynamic  signal  analyzer,  step  by  step  details  are  given  for  setup 
and  use  of  the  dynamic  signal  analyzer  for  spectral  measurements  in 
Appendix  C. 
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With  settings  on  the  signal  conditioning  units  adjusted  as  shown  in  Table 
6,  frequency  spectra  were  obtained  between  0.1  and  25  Hz  and  between  0.1 
and  10000  Hz. 

2.  Probe  Calibration  and  Measurement  Procedures 

The  probe  was  calibrated  in  the  test  section  of  the  wind  tunnel 
located  in  the  Department  of  Mechanical  Engineering  laboratory  of  the 
Naval  Postgraduate  School.  From  the  calibration,  N=0.45,  B=3. 750916,  and 
EOC=2. 653925,  where  the  equation  given  by: 


E2-EOC 2 
_ 


\l/n 

/ 


is  valid  when  the  sensor  is  subjected  to  forced  convection.  For  Red  <0.07 
mixed  convection  affects  sensor  response  and  the  relationship  between 
velocity  and  voltage  is  given  by  the  equation: 

Y=A0+A1  *X+A2  *X2+A3  *X3 

where  Y=UN,  X=E2-EOC2,  A0=-0. 68058,  A1=0. 99448,  A2=-0. 24088  and 

A3-0. 024986. 

With  the  experimental  arrangement,  longitudinal  velocity 
fluctuation  magnitudes  are  determined  from  digital  processing  of  hot-wire 
bridge  voltage  signals  and  from  signals  read  from  the  DISA  56N25  RMS  unit. 
For  the  former,  20000  voltage  samples  were  obtained.  For  the  latter,  the 
normalized  streamwise  velocity  fluctuation  magnitudes  are  calculated  using 
the  equation  from  given  by: 

o',  2  j>' 

u  "  n  fa2 -A) 


[Ref.  17]. 
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APPENDIX  B  -  SOFTWARE  DIRECTORY 


Programs  written  in  Hewlett  Packard  BASIC  4.0  language  for  the  HP  300 
Series  computers  utilized  in  these  studies  are  listed.  A  description  of 
wha>u.  the  program  does,  significant  user  inputs  required  and  significant 
output  is  provided. 


A.  GENERAL  PROGRAMS 

DEAN15M  Calculates  the  Dean  number  as  described  in  Chapter  II. 
user  input 

•  barometric  pressure 

•  ambient  temperature 

•  orifice  pressure  drop 

output 

•  expansion  coefficient  Y 

•  flow  coefficient  K 

•  mass  flow  rate 

•  channel  Reynolds  number 

•  Dean  number 

•  bulk  velocity 


B.  FIVE-HOLE  PRESSURE  PROBE  PROGRAMS 

C4  Data  file  of  calibration  coefficients:  pitch,  yaw,  Cpy,  C^, 

C  ..  C  .  and  C  .  . 
pt'  ps'  pts 

£4  Data  file  of  calibration  pressures:  pitch,  yaw,  P^,  P2,  P3,  P4, 
and  P5. 

FIVEHOLE  Acquires  pressure  measurements  at  each  of  the  five  ports 
and  directs  probe  positioning. 


user  input 

•  date/time 

•  estimated  dean  number 

•  orifice  size 

•  number  of  spanwise  and  radial  points 

•  spanwise  and  radial  resolutions 

•  spanwise  and  radial  starting  positions 

•  wait  time  for  pressure  to  stabilize 

•  number  of  transducers 

•  number  of  samples  per  channel 
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and  curvature  of  115°  on  the  right  side  of  the  channel.  Still  photography 
timed  sequences  showing  a  collapsing  occurrence  and  a  splitting  occurrence 
are  presented  also.  For  each  De,  these  data  were  obtained  from 
spanwise/radial  plane  observations  which  were  video  tape  recorded  at 
angles  from  85°  to  135°  from  the  start  of  curvature. 

A.  SURVEYS  OF  TIME-AVERAGED  QUANTITIES 

Surveys  of  time-averaged  quantities  determined  from  five-hole  pressure 
probe  measurements  are  presented  in  figures  17  through  34.  Figures  17 
through  19  present  total  pressure  contours.  Figures  20  through  22  present 
streamwise  velocity  contours.  Figures  23  through  25  present  streamwise 
velocity  perturbations.  Figures  26  through  28  present  streamwise 
vorticity  components.  Figures  29  through  31  present  radial  vorticity 
components.  Figures  32  through  34  present  spanwise  vorticity  components. 

1.  Total  Pressure 

In  Figures  17  through  19,  surveys  of  Pamb-Ptotal  (in  Pascals) 
representing  ambient  pressure  minus  measured  total  pressure  are  given. 
The  higher  the  index  in  the  flow,  the  higher  the  difference  in  pressure, 
i.e.,  the  higher  the  vacuum  measured.  Pressure  deficits  near  the  concave 
surface  in  the  form  of  humps  are  upwash  regions  between  counter-rotating 
vortex  pairs.  At  De=75  the  flow  shows  evidence  of  portions  of  at  least 
five  deficit  regions.  At  De=125.1  the  number  of  deficits  is  reduced  to 
two.  The  highest  deficits  occur  at  De=175.2.  At  De>300,  the  surveys  show 
reasonable  spanwise  uniformity. 

2 .  Streamwise  Velocity 

Surveys  of  streamwise  velocity  normalized  by  bulk  velocity  are 
displayed  in  Figures  20  through  22.  These  show  velocity  deficits  near  the 
concave  wall  which  are  nearly  identical  qualitatively  to  the  total 
pressure  deficits  in  Figures  17  through  19.  These  deficits  are  less 
intense  than  ones  measured  in  the  curved  portion  of  the  channel  by  Fields 
(Ref.  15]  because  the  mechanism  to  produce  the  vortices  and  accompanying 
secondary  flows  (centrifugal  instabilities)  is  not  present. 


19 


3 .  Velocity  Perturbation 


In  order  to  show  streamwise  velocity  perturbations,  curved  channel 
Poiseuille  flow  is  subtracted  from  streamwise  velocity.  Such 
perturbations  are  presented  in  Figures  23  to  25.  At  De=100.3  a  regular 
perturbation  pattern  is  present  across  the  span  of  the  measurement  plane. 
The  perturbations  become  spanwise  uniform  for  De>300. 

4.  Vorticity  Components 

Vorticity  components  in  Figures  26  through  34  have  positive  and 
negative  signs.  Positive  vorticity  is  contained  within  the  solid  contour 
lines  while  negative  vorticity  is  contained  within  the  dashed  contour 
lines.  As  the  plots  are  arranged  so  that  streamwise  flow  moves  into  the 
page,  positive  streamwise  vorticity  is  clockwise. 

a.  Streamwise  Vorticity 

Streamwise  vorticity  contours  are  shown  in  Figures  26  through 
28.  The  streamwise  vorticity  distributions  correlate  nicely  to  total 
pressure  and  streamwise  velocity  such  that  velocity/pressure  deficits  are 
present  between  the  two  vortices  within  each  counter-rotating  pair. 
Counter-rotating  vortex  pairs  are  clearly  evident  at  the  locations  of  the 
two  largest  total  pressure  deficits  for  De=100.3.  At  De=175.2  two  pairs 
of  well  defined  vortices  again  correspond  to  tie  locations  of  the  velocity 
deficits.  At  De  greater  than  about  300  streamwise  vorticity  contours  tend 
toward  spanwise  uniformity. 

b.  Radial  Vorticity 

Radial  vorticity  contours  are  shown  in  Figures  29  through  31. 
As  De  increases,  the  flow  becomes  spanwise  periodic  at  De=75,  until  De>300 
when  the  flow  becomes  spanwise  uniform. 

c.  Spanwise  Vorticity 

Spanwise  vorticity  contours  are  shown  in  Figures  32  through 
34.  These  contours  are  consistent  with  the  total  pressure  contours  and 
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contours  of  streamwise  velocity,  where  negative  vorticity  components 
closely  match  the  velocity/pressure  deficits.  Flow  is  spanwise  uniform  at 
De=50.4  and  for  De>300. 

B.  TIME-AVERAGED  AND  SPATIALLY  RESOLVED  FLOW 

PROPERTIES  NEAR  THE  CONVEX  SURFACE 

Results  presented  in  this  section  were  obtained  using  an  array  of  hot- 
film  gages  mounted  on  the  convex  surface  of  the  channel  (y/d=1.0)  120° 
from  the  start  of  curvature.  Figure  35  shows  the  locations  of  the  12  hot- 
film  elements  with  respect  to  vortex  pairs.  Table  1  gives  probe  locations 
with  respect  to  individual  vortices.  One  expects  results  from  probes  near 
similar  vortex  characteristics  to  be  correlated  with  each  other.  Figures 
36  through  43  give  correlations  versus  z/d  referenced  to  probe  2  for 
probes  2,  4,  6,  8,  10  and  12,  for  De  from  50.25  to  399.8.  Figure  44  is 
then  a  3-D  composite  of  Figures  36  through  43.  Figures  45  through  52  give 
correlations  versus  z/d  referenced  to  probe  12  for  probes  2,  4,  6,  8,  10 
and  12,  for  De  from  50.25  to  399.8.  Figure  53  is  then  a  3-D  composite  of 
Figures  45  through  52. 

Referring  to  Figures  44  and  53  shows  that,  in  all  cases,  correlation 
magnitudes  decrease  as  the  distance  from  the  reference  probe  increases. 
Reference  probes  are  probe  2,  located  at  z/d=4.73,  for  Figure  44  and  probe 
12,  located  at  z/d=7.23  for  Figure  53.  In  Figure  44,  sign  variations  of 
correlations  versus  z/d  are  evident  at  De  from  50.25  to  299.6.  These 
variations  give  evidence  of  regions  of  positive  and  negative  correlation 
between  different  parts  of  vortex  pair  structures.  For  De=100,  probe  2  is 
located  directly  over  a  downwash  region.  Probe  4,  the  nearest  active 
probe,  is  located  directly  over  an  upwash  region  and  its  signal  correlates 
negatively  with  the  signal  from  probe  2.  At  De>299.6,  correlations  have 
fully  turbulent  characteristics,  dropping  to  values  which  do  not  vary 
significantly  with  z/d. 

Similar  qualitative  trends  are  evident  in  Figure  53.  For  De=100, 
probe  12  is  located  over  a  region  of  negative  vorticity.  Probe  10,  the 
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nearest  active  probe,  is  located  over  a  region  of  positive  vorticity  and 
its  signal  correlates  negatively  with  the  signal  from  probe  12.  At 
De>250.4,  correlations  have  fully  turbulent  characteristics,  dropping  to 
values  which  do  not  vary  significantly  with  z/d. 

Figure  54  shows  the  average  (over  all  gages  employed)  normalized  rms 
voltage  versus  De  from  50.25  to  399.8.  This  curve  indicates  an  increase 
in  Vrms/Vmean  as  De  increases  from  150  to  200.  At  De  higher  than  300,  the 
variations  with  De  have  lower  slopes. 

Figure  55  shows  scaled,  normalized  rms  voltage  versus  z/d  for  each  De 
surveyed.  The  normalized  rms  voltage  for  each  channel  is  scaled  by  the 
normalized  rms  voltage  measured  at  De=399.8.  Thus  the  scaled  values  at 
this  De  are  equal  to  1.0.  These  curves  show  that  spatially  varying 
structures  occur  at  De=149.9,  200.7,  250.4,  299.6  and  350.5. 

Length  scales  shown  in  Figure  56  represent  integrations  of  the 
correlations  with  respect  to  z/d.  These  evidence  spatially  varying 
structures  between  De=100  and  De=250. 

The  previous  results  were  discussed  for  an  active  probe  spacing  of 
6.35  mm.  For  a  second  set  of  measurements,  active  probe  spacing  was 
reduced  to  3.18  mm.  Results  are  now  presented  which  were  obtained  UBing 
the  array  of  hot-film  gages  with  the  smaller  probe  spacing.  Figures  57 
through  64  give  correlations  versus  z/d  referenced  to  probe  1  for  probes 
1,  2,  3,  4,  5  and  6,  for  De  from  49.2  to  400.3.  Figure  65  is  then  a  3-D 
composite  of  Figures  57  through  64.  Figures  66  through  73  give 
correlations  versus  z/d  referenced  to  probe  6  for  probes  1,  2,  3,  4,  5  and 
6,  for  De  from  49.2  to  400.3.  Figure  74  is  then  a  3-D  composite  of 
Figures  66  through  73. 

Referring  to  Figures  65  and  74  shows  that,  in  all  cases,  correlation 
magnitudes  decrease  as  the  distance  from  the  reference  probe  increases. 
Reference  probes  are  probe  1,  located  at  z/d=4.48,  for  Figure  65  and  probe 
6,  located  at  z/d=5.73  for  Figure  74.  In  Figure  65,  sign  variations  of 
correlations  versus  z/d  are  evident  at  De  from  49.2  to  200.9.  These 
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variations  give  evidence  of  regions  of  positive  and  negative  correlation 
between  different  parts  of  vortex  pair  structures.  For  De=100,  probe  1  is 
located  directly  over  a  region  of  positive  vorticity.  Probe  2,  the 
nearest  active  probe,  is  located  directly  over  a  downwash  region  and  its 
signal  correlates  negatively  with  the  signal  from  probe  1.  At  Dea249.6, 
correlations  have  fully  turbulent  characteristics,  dropping  to  values 
which  do  not  vary  significantly  with  z/d. 

Similar  qualitative  trends  are  evident  in  Figure  74.  For  De=100, 
probe  6  is  located  over  a  downwash  region.  Probe  5,  the  nearest  active 
probe,  is  located  over  a  region  of  positive  vorticity  and  its  signal 
correlates  negatively  with  the  signal  from  probe  6.  At  Dea350.2, 
correlations  have  fully  turbulent  characteristics,  dropping  to  values 
which  do  not  vary  significantly  with  z/d. 

Figure  75  shows  the  average  (over  all  gages  employed)  normalized  rms 
voltage  versus  De  from  49.2  to  400.3.  This  curve  indicates  an  increase  in 
Vrms/Vmean  as  De  increases  from  200  to  250.  At  De  higher  than  250,  the 
variations  with  De  have  lower  slopes. 

Figure  76  shows  scaled,  normalized  rms  voltage  versus  z/d  for  each  De 
surveyed.  The  normalized  rms  voltage  for  each  channel  is  scaled  by  the 
normalized  rms  voltage  measured  at  De=400.3.  Thus,  as  for  Figure  55,  the 
scaled  values  at  this  De  are  equal  to  1.0.  These  curves  show  that  the 
spatially  varying  structures  occur  at  De=100.1,  150.7,  200.9,  249.6,  299.2 
and  350.2. 

Length  scales  shown  in  Figure  77  represent  integrations  of 
correlations  with  respect  to  z/d.  These  evidence  spatially  varying 
structures  between  De=100  and  De=200. 

C.  SPECTRAL  AND  LONGITUDINAL  VELOCITY  FLUCTUATION 

VARIATIONS  WITH  DEAN  NUMBER 

Figure  16  shows  the  five  hot-wire  probe  locations  where  surveys  were 
taken  in  the  spanwise/radial  plane  120°  from  the  start  of  curvature.  At 
each  of  these  locations,  time-averaged  streamwise  velocity  and 
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longitudinal  turbulence  intensity  were  measured  using  a  single  hot-wire 
probe  at  De  from  50  to  450.  The  five  locations  were  chosen  to  obtain 
information  on  different  parts  of  the  vortex  pair  structure.  As  shown  in 
Figure  16,  position  1  is  located  near  the  concave  wall  at  the  base  of  an 
upwash  region  (with  respect  to  the  concave  wall).  Position  2  is  located 
in  an  upwash  region  at  the  channel  radial  centerline.  Position  3  is 
located  near  the  concave  wall  near  the  base  of  a  downwash  region  (with 
respect  to  the  concave  wall).  Position  4  is  located  near  the  concave  wall 
directly  under  a  vortex.  Position  5  is  located  in  a  downwash  region  at 
the  channel  radial  centerline. 

The  streamwise  mean  velocity  and  longitudinal  turbulence  intensity 
were  each  determined  two  ways.  With  a  digital  method,  the  instantaneous 
streamwise  velocity  was  calculated  for  each  of  20000  samples.  The  average 
of  the  20000  instantaneous  velocities  is  then  the  mean  streamwise 
velocity.  The  longitudinal  velocity  fluctuation  intensity  was  then 
determined  as  the  root-mean-square  of  the  instantaneous  longitudinal 
velocities.  With  an  analog  method,  mean  velocity  and  longitudinal 
turbulence  intensity  were  calculated  based  on  mean  and  rms  voltages 
integrated  over  10  second  time  intervals. 

Figures  78  and  79  show  plots  of  the  average  streamwise  velocity  using 
digital  and  analog  techniques,  respectively,  versus  De  for  the  five  probe 
positions  listed  in  Table  5.  The  digital  and  analog  streamwise  velocities 
coincide  nicely.  All  positions  show  a  velocity  increase  that  is  generally 
linear  with  De.  The  behavior  of  the  velocity  for  positions  1,  3  and  4 
varies  because  spatial  locations  of  upwash  regions  vary  with  De.  The 
proximity  of  positions  1,  3  and  4  to  the  concave  wall  causes  associated 
velocity  magnitudes  to  be  generally  lower  than  for  positions  2  and  5. 

Figures  80  and  81  show  longitudinal  turbulence  intensities  determined 
using  digital  and  analog  techniques,  respectively,  versus  De  for  the  five 
probe  positions.  The  digital  and  analog  data  coincide  with  each  other 
reasonably  well.  For  De<160,  the  intensities  are  very  small  for  all  five 
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probe  positions  .  At  De  higher  than  160,  intensities  generally  increase 
with  De  approximately  linearly.  Positions  2  and  5  have  lower  intensities 
because  they  are  farther  away  from  the  effects  of  the  concave  surface  than 
positions  1,  3  and  4. 

Normalized  longitudinal  turbulence  intensities  are  shown  in  Figures  82 
and  83.  Digital  and  analog  data  coincide  well  at  most  all  De,  except  for 
values  of  De=50  when  the  digitally  obtained  values  are  incorrectly 
indicated  to  be  too  high.  For  De  higher  than  160,  longitudinal  turbulence 
intensity  increases  with  De.  The  variations  at  position  4  are  due  to 
time-averaged  spatial  location  variance  of  the  upwash  regions  with  De. 
Over  this  range  of  De,  values  for  probe  positions  2  and  5  are  lower  than 
others  because  positions  2  and  5  are  farther  away  from  the  concave  wall. 

Figure  84  gives  settings  used  on  the  dynamic  signal  analyzer  to  obtain 
power  spectra  over  frequencies  from  0.1  to  25  Hz.  Summaries  of  observed 
local  frequency  maxima  are  given  in  Figures  85  through  89  for  De=50  to 
De=200  for  each  of  the  five  probe  positions  1  through  5,  respectively. 
Frequencies  of  1  Hz  to  2  Hz,  5  Hz  to  6  Hz,  8  Hz,  and  20  Hz  to  21  Hz  exist 
for  all  five  probe  positions.  An  additional  frequency  of  about  3  Hz 
exists  for  probe  positions  3,  4  and  5.  Figures  90  and  91  are  two  examples 
of  power  spectra  used  to  obtain  such  information.  Figures  85  through  89 
were  compiled  from  data  presented  in  Figures  92  through  151,  which  show 
individual  frequency  spectra  for  all  five  probe  positions  at  De  from  50  to 
450. 

Figure  152  gives  settings  used  on  the  dynamic  signal  analyzer  to 
obtain  high  frequency  power  spectra  over  frequencies  from  0.1  Hz  to  10.0 
kHz.  Figure  153  shows  the  average  of  five  power  spectra  at  probe  position 
2  at  De=170.  Clearly  evident  is  a  principal  peak  between  90  Hz  and  100 
Hz.  Also  apparent  are  harmonic  peaks  at  about  180  Hz,  280  Hz  and  400  Hz. 
Such  variations  are  believed  to  be  due  to  twisting  [Ref.  18].  Figure  154 
shows  a  plot  of  principal  peak  magnitudes  versus  De  for  all  five  probe 
positions  at  De  from  150  to  200.  In  this  figure,  peak  magnitudes  are 


shown  to  be  dependent  on  De  and  probe  position,  where  peak  magnitudes 
generally  increase  with  increasing  De.  The  highest  peak  magnitudes 
observed  are  for  De=180  at  probe  position  4  and  for  De=200  at  probe 
position  1.  Such  behavior  suggests  higher  instability  near  the  upwash 
region  close  to  the  concave  surface.  Figure  154  was  compiled  from  data 
shown  in  Figures  155  through  214  which  are  individual  power  spectra  for 
all  five  probe  positions  at  De  from  50  to  450. 

D.  ANALYSIS  OF  APPEARANCE  AND  DISAPPEARANCE  PHENOMENA 

FROM  FLOW  VISUALIZATION  RESULTS 

Appearances  and  disappearances  of  already  developed  vortex  pairs  along 
with  accompanying  alterations  of  spanwise  vortex  pair  spacing  are 
discussed  by  Ligrani  et  al .  [Ref.  6]  At  the  time  of  this  writing,  there 
is  evidence  of  two  types  of  appearance  events  and  four  types  of 
disappearance  events.  Appearance  events  are  categorized  as  growth  or 
splitting.  Growth  events  occur  when  vortex  pairs  emerge  from  the  flow 
region  near  the  concave  wall.  Splitting  events  occur  when  one  vortex  pair 
is  rearranged  into  two  vortex  pairs.  Splitting  is  also  observed  after  a 
vortex  pair  decreases  in  size  with  simultaneous  emergence  of  two  new 
vortex  pairs  nearby  from  flow  near  the  concave  Burface. 

Disappearance  events  are  categorized  as  engulfment,  cancellation, 
merging  and  collapse.  Engulfment  refers  to  the  absorbtion  of  a  vortex 
pair  by  another  vortex  pair  where  contact  occurs  between  vortices  with  the 
same  sign  of  vorticity.  Cancellation  occurs  when  vortices  of  tho  opposite 
sign  collide.  Typically,  these  vortex  pairs  rock  slightly  toward  each 
other  as  they  move  together.  With  merging,  a  smaller  vortex  pair  merges 
into  the  upwash  region  of  a  larger  pair.  Collapse  events  are  evidenced  by 
the  disappearance  of  a  pair  into  the  flow  region  near  the  concave  surface. 

Figure  215  shows  a  new  type  of  splitting  event  which  was  unrecognized 
prior  to  the  present  work.  Here,  the  center  vortex  pair  in  the  field  of 
view  shrinks  in  size  until  it  disappears  into  the  flow  near  the  concave 
surface  while  two  new  vortex  pairs  emerge  simultaneously  from  the  base  of 
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the  shrinking  vortex  pair.  The  sequence  starts  at  the  top  with  frame  1. 
Each  photo  shows  the  full  channel  height,  12.7  mm,  and  a  channel  width  of 
50.8  mm,  and  gives  a  spanwise/radial  plane  view  with  the  concave  surface 
at  the  bottom  and  streamwise  flow  into  the  plane  of  the  page.  In  frames 

I  through  5  the  center  vortex  pair  begins  to  shrink.  In  frame  6,  dark 
spots  begin  to  form  on  each  side  at  the  base  of  the  shrinking  vortex  pair. 
Both  spots  rapidly  grow  in  frames  7  and  8  while  the  original  vortex  pair 
shrinks  even  more.  In  frame  9  the  original  vortex  pair  is  no  longer 
evident.  In  frame  10  the  two  new  vortex  pairs  are  fully  developed.  In 
frames  11  and  12  the  larger  of  the  two  new  vortex  pairs  rocks  to  the  left. 
Also  in  frames  11  and  12,  the  new  vortex  pair  on  the  left  begins  to  shrink 
toward  the  concave  wall.  The  larger  of  the  two  new  vortex  pairs  is  now 
displaced  to  the  left. 

Figure  216  shows  a  collapsing  type  of  disappearance  event.  Here,  the 
second  vortex  pair  from  the  right  will  shrink  until  it  disappears  into  the 
flow  near  the  concave  wall.  Frames  1  through  4  show  the  smallest  vortex 
pair  present  growing  to  full  size.  In  frames  5  and  6,  it  twists  slightly 
counterclockwise  and  begins  to  shrink.  Frames  7  through  10  then  show  the 
rapid  collapse  into  the  flow  near  the  concave  wall  by  the  small  vortex 
pair  and  then  the  translation  to  the  left,  as  it  fills  the  void  left 
behind,  by  the  large  vortex  pair  at  the  right  side  of  the  frames.  Frames 

II  and  12  also  show  the  center  vortex  pair  enlarging  to  occupy  a  larger 
space. 

Figure  217  shows  a  time  history  of  appearance  and  disappearance  events 
observed  over  an  80  second  interval  at  De=100  and  115°  from  the  start  of 
curvature  as  observed  on  the  right  side  of  the  channel.  Figure  218  gives 
magnified  views  of  three  regions  of  Figure  217.  The  scale  on  these  three 
time  lines  allows  determination  of  the  time  intervals  between  discrete 
events.  The  time  lines  show  that  most  appearance  events  are  followed  by 
disappearance  events,  and  vice-versa. 
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Frequencies  of  appearance  and  disappearance  events  are  given  in  Tables 
2,  3  and  4  in  dimensional  and  non-dimensional  forms.  The  tabular  data 
from  Tables  3  and  4  is  also  shown  in  Figures  219  to  232.  Figures  219,  221 
and  223  present  frequencies  of  appearance  and  disappearance  events  versus 
De  for  angles  of  85®,  95°  and  105°  from  the  start  of  curvature 
respectively.  Figures  220,  222  and  224  present  the  same  information  in 
non-dimensional  form  as  in  Figures  219,  221  and  223  respectively.  Figures 
225  and  226  present  frequencies  of  appearance  and  disappearance  events 
versus  0  at  De=75  and  De=100  respectively.  Figure  227  shows  frequencies 
of  total  appearance  and  total  disappearance  events  versus  De  for  angles 
from  85®  to  135®  after  the  start  of  curvature.  Figures  228  and  229 
provide  curve  identification  for  Figure  227.  Figures  230,  231  and  232 
present  the  same  information  in  non-dimensional  form  as  in  Figures  227, 
228  and  229  respectively.  Frequencies  of  appearance  and  disappearance 
events  generally  increase  with  increasing  De  for  Figures  219  through  224 
and  227  through  232.  Figure  225  shows  that  frequencies  of  appearance  and 
disappearance  events  generally  increase  with  increasing  0  at  De=75. 
Figure  226  shows  that  frequencies  of  appearance  and  disappearance  eventB 
generally  decrease  with  increasing  0  at  De=100.  In  all  cases,  growth  type 
appearance  events  and  collapse  type  disappearance  events  occur  most  often. 
In  Figures  219  and  225,  local  minima  of  total  appearance  and  total 
disappearance  event  frequencies  occur  at  De=125  and  85°  after  the  start  of 
curvature  and  De=75  and  115®  after  the  start  of  curvature  respectively. 
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V.  SUMMARY  AND  CONCLUSIONS 

Instabilities  in  curved  channel  flow  were  studied  over  conditions 
ranging  from  fully  laminar  to  fully  turbulent.  Primary  instabilities  in 
the  form  of  regularly  spaced  Dean  vortices  develop  in  the  flow  due  to 
centrifugal  effects.  After  the  primary  instabilities,  secondary 
instabilities  develop  which  affect  the  behavior  of  the  vortices.  These 
secondary  instabilities  take  the  form  of  undulation  and  twisting  of  the 
vortices,  plus  phenomena  where  vortex  pairs  appear  and  disappear.  Two 
curved  channels  with  40  to  1  aspect  ratios  were  used  to  obtain  information 
on  these  instabilities. 

Time-averaged  total  pressure,  mean  velocity,  velocity  perturbation  and 
vorticity  component  surveys  were  obtained  using  a  miniature  five-hole 
pressure  probe  at  De  from  50.4  to  424.8  at  a  streamwise  location  in  the 
second  straight  section  located  downstream  of  the  curved  portion  of  the 
channel.  At  De<75,  results  show  that  flow  is  spanwise  uniform.  As  De 
increases,  results  provide  evidence  of  vortex  pairs  located  at  spanwise 
locations  across  the  measurement  plane.  At  De  higher  than  about  300,  flow 
again  becomes  spanwise  uniform. 

Time-averaged  and  spatially  resolved  flow  correlations  were  obtained 
using  a  12  element  hot-film  probe  located  120°  from  the  start  of  curvature 
at  y/d=1.0  (flush  with  the  convex  wall)  for  De  from  49.2  to  400.3.  Six 
probe  elements  were  used  in  two  different  probe  configurations  with 
element  spacings  of  6.35  mm  and  3.18  mm  respectively.  In  all  cases, 
correlation  magnitudes  decrease  as  the  spanwise  distance  from  the 
reference  probe  increases.  Sign  variations  of  correlations  versus  z/d 
evidence  regions  of  positive  and  negative  correlation  between  different 
parts  of  vortex  pair  structures.  For  the  first  configuration,  average 
(over  all  gages  employed)  normalized  rms  voltage  increases  as  De  increases 
from  De=150  to  De=200.  Scaled,  normalized  local  rms  voltages  versus  z/d 
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show  spatially  varying  structures  at  De  from  150  to  350.  Length  scales 
indicate  spatially  varying  structures  at  De  from  100  to  250. 
Qualitatively,  the  second  probe  configuration  is  in  close  agreement  with 
the  first  configuration. 

Power  spectra,  streamwise  velocity  and  longitudinal  turbulence 
intensity  surveys  were  obtained  using  a  hot-wire  probe  at  De  from  50  to 
450.  Each  survey  was  obtained  at  each  of  five  positions  in  the 
spanwise/radial  plane  120°  from  the  start  of  curvature.  Digital  and 
analog  methods  were  used  to  determine  streamwise  velocity  and  longitudinal 
turbulence  intensity.  Results  from  both  methods  are  in  close  agreement. 
Streamwise  velocity  generally  increases  linearly  over  the  range  of  De. 
The  behavior  of  velocity  for  positions  near  the  concave  wall  varies 
because  spatial  locations  of  upwash  regions  with  respect  to  the  concave 
wall  vary  with  De.  Longitudinal  turbulence  intensity  increases  with  De 
higher  than  160.  Positions  near  the  concave  wall  have  higher  intensities 
than  the  positions  at  the  radial  centerline.  Observed  local  frequency 
maxima  from  low  frequency  power  spectra  evidence  frequencies  of  1  Hz  to  2 
Hz,  5  Hz  to  6  Hz,  8  Hz,  and  20  Hz  to  21  Hz  for  all  five  probe  positions. 
Local  frequency  maxima  also  exist  at  3  Hz  for  probe  positions  away  from 
upwash  regions.  For  high  frequency  power  spectra,  principle  peak 
magnitudes  increase  with  increasing  De.  Highest  peak  magnitudes  occur  at 
positions  near  the  base  of  the  upwash  region  which  evidence  twisting  [Ref. 
18] . 

Frequencies  of  events  where  vortex  pairs  appear  and  disappear  are 
obtained  for  De  from  60  to  150.  Time  lines  showing  these  events  are  shown 
at  De=100.  Still  photography  time  lapse  sequences  are  presented  showing 
the  collapse  of  a  vortex  pair  and  one  type  of  splitting  of  a  vortex  pair. 
These  data  were  obtained  for  each  De  from  spanwise/radial  plane 
observations  which  were  video  tape  recorded  at  angles  from  85°  to  135° 
from  the  start  of  curvature.  Time  line  results  indicate  that  should  an 
appearance  event  occur,  the  likelihood  is  that  a  disappearance  event  will 


occur  next,  and  vice-versa.  Frequencies  of  total  appearance  and  total 
disappearance  events  are  nearly  equal  in  all  cases.  Frequencies  generally 
increase  with  increasing  De.  Growth  type  of  appearance  events  and 
collapse  type  of  disappearance  events  occur  most  often.  Highest  total 
appearance  and  total  disappearance  frequencies  observed  are  at  De=125  and 
95°  after  the  start  of  curvature.  Local  minima  of  total  appearance  and 
total  disappearance  event  frequencies  occur  at  De=125  and  85°  after  the 
start  of  curvature  and  at  De=75  and  115°  after  the  start  of  curvature. 
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APPENDIX  A  -  UNCERTAINTY  ANALYSIS 


The  uncertainty  analysis  presented  below  for  sections  A.  and  B.  is 
taken  from  Baun  [Ref.  9].  The  present  study  utilizes  much  of  the  same 
peripheral  equipment  and  studies  flow  situations  similar  to  those  studied 
by  Baun  in  the  same  channels. 

A.  DEAN  NUMBER  UNCERTAINTY 

The  Dean  number  is  a  function  of  orifice  pressure  drop.  The 
uncertainty  for  the  pressure  drop  across  the  orifice  is  analyzed  at 
De=150. 

1.  Pressure  Drop  Uncertainty 
6APO(.=±0.001  in.  HjO 

2.  Mass  Flow  Rate  Uncertainty 
5AQr=±0 . 012  in2 

5K=±0 . 001 
5Y=±0.02 

6p=±0. 002  lbm/ft3 

6m=±0.0004  lbm/sec  (95%  confidence  level) 

6Ach=±0.0052  in2 

6De=±4.4  (95%  confidence  level) 


B.  FIVE-BOLE  PROBE  VELOCITY  MEASUREMENT  UNCERTAINTY 

Five-hole  probe  velocity  component  uncertainties  are  derived  from 
calibration  uncertainties  and  measurement  uncertainties.  Values  at 
De=100.3,  y/d=0.5  and  z/d=0.0  are  listed  below.  The  probe  calibration 
uncertainty  is  for  yaw=-10°  and  pitch=0.0°. 

1.  Validyne  Transducer  Calibration  Uncertainty 

5Pnanon*ter=±0-0°l  in’  H2° 

6E=i0.0007  Volts 
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SC=±0. 000072  in.  HgO/Volt  (bias  error  is  disregarded  because 

pressure  differentials  are  calculated) 

7.  Probe  Calibration  Uncertainty 

5E=±0.0004  Volts 

6P.|=±0. 0006  in.  HjO 

6P2=±0.0001  in.  H20 

5P3=±0.0005  in.  H20 

5P4=±0.0002  in.  H20 

5P5=±0.0003  in.  H20 

6PS=±0. 00004  in.  H20 

6Pt=±0.0007  in.  HgO 

6P=±0. 00015  in.  H20 

5Cpy=±0.025 

6Cpp=±0.  Oil 

6C  f  =*0. 045 
pts 

6Ky=±0.004  degrees'1 
6Kp=±0.004  degrees'1 
6Kts=±0.003  degrees  _1 

3.  Curved  Channel  Measurement  Uncertainty 

5E=±0 . 0007  Volts 
6C=±0. 000072  in.  H20/Volt 
6P^=±0. 0002  in.  H20 
5P2=±0.0002  in.  H20 
6P3=±0.0002  in.  H20 
6P4=±0.0002  in.  H20 
5P5=±0.0002  in.  HzO 
5P=±0.0001  in.  H20 
5Cpy=±0.134 
6cpp=±0. 14 1 
6yaw=±l. 33° 

6pitch=±1.4° 
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5 (Pt-Ps) =±0.0004  in.  HgO 
6U=±0.063  m/s 
6Ug=±0.062  m/s 
6Ur=±0.062  m/s 
5UZ=±0.030  m/s 

C.  HOT-WIRE  PROBE  VELOCITY  MEASUREMENT  UNCERTAINTY 

The  uncertainty  associated  with  the  hot-wire  probe  velocity  is  derived 
from  probe  calibration  uncertainty  and  the  uncertainty  of  voltages 
produced  by  the  probe  and  bridge  network. 

1.  Probe  Calibration  Uncertainty 
6B=±0. 0216  Volts2  s-45  m'-;e 
6EQ=±0. 0005  Volts 

6EOC=±0 . 0005  Volts 
6EOM=±0.0005  Volts 
6AP=±0.0005  in.  H20 

2.  Curved  Channel  Measurement  Uncertainty 
6Eo=±O.0005  Veits 

6UX=±0.048  m/s 
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APPENDIX  B  -  SOFTWARE  DIRECTORY 


Programs  written  in  Hewlett  Packard  BASIC  4.0  language  for  the  HP  300 
Series  computers  utilized  in  these  studies  are  listed.  A  description  of 
what  the  program  does,  significant  user  inputs  required  and  significant 
output  is  provided. 


A.  GENERAL  PROGRAMS 

DEAN15M  Calculates  the  Dean  number  as  described  in  Chapter  II. 


user  Inoat 

•  barometric  pressure 

•  ambient  temperature 

•  orifice  pressure  drop 

output 

•  expansion  coefficient  Y 

•  flow  coefficient  K 

•  mass  flow  rate 

•  channel  Reynolds  number 

•  Dean  number 

•  bulk  velocity 


B.  FIVE-HOLE  PRESSURE  PROBE  PROGRAMS 


C4 


Data  file  of  calibration  coefficients:  pitch,  yaw,  C  ,  C  , 

r  /  rr 


Cpt'  Cps'  and  Cpts* 

P4  Data  file  of  calibration  pressures:  pitch,  yaw,  P1 ,  P2,  P3,  P^, 


and  P5. 

FIVEHOLE  Acquires  pressure  measurements  at  each  of  the  five  ports 
and  directs  probe  positioning. 
user  input 


•  date/time 

•  estimated  dean  number 

•  orifice  size 

•  number  of  spanwise  and  radial  points 

•  spanwise  and  radial  resolutions 

•  spanwise  and  radial  starting  positions 

•  wait  time  for  pressure  to  stabilize 

•  number  of  transducers 

•  number  of  samples  per  channel 
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•  number  of  standard  deviations 

•  calibration  coefficients 

•  ambient  pressure 

output 

•  FIVxxx  data  file  with:  y,  z,  C  ,  C  ,  P.,  P,  De 

•  FIVPxxx  data  file  with  y,  z,  P^7  P2»  P3»  Pg,  p*  De 

•  average  De  : 

PAD JUST  Corrects  spatial  resolution  described  in  Chapter  III. 

user  input 

•  FIVPxxx  data  file 

•  yien 

•  zlen 
output 

•  FIVPxxxA  data  file  with  y,  z,  Cy,  Cp,  P1f  P,  De 

ORIENT  Orients  the  five-hole  probe  in  the  flow  such  that  P2=P3  at 
yaw=0°  and  pitch=0° 
user  input 

•  number  of  samples  per  channel 

•  number  of  standard  deviations 

•  transducer  calibration  coefficients 

•  barometric  pressure 

output 


•  mean  velocity 

•  C 

• 

•  p;,  p2,  p3,  p 


4'  P5 


VELOCITY 1 

user  input 


Computes  the  three  mean  velocity  components 


•  C4  data  file 

•  P4  data  file 

•  number  of  yaw  and  pitch  points 

•  FIVPxxxA  data  file 

•  pitch  angle  offset  (4.3°) 

•  yaw  angle  offset  (1.0°) 

•  y  and  z  direction  downwash  velocity  corrections  (0.011,  0.011) 
output 

•  Vxxx  data  file  with  y,  z,  P{,  U,  Ue,  Ur,  Uz 
Uxl  Plots  mean  streamwise  velocity  contours 

user  input 

•  Vxxx  data  file 
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•  number  and  spacing  of  contour  levels 
output 

•  contour  plot  with  level  table 

PT0T1  Plots  Pamb  minus  Ptotal  contours 
user  input 

•  Vxxx  data  file 

•  number  and  spacing  of  contour  levels 
output 

•  contour  plot  with  level  table 
CCPF1  Calculates  CCPF 

user  input 

•  Vxxx  data  file 
output 

•  PFxxx  data  file  with  y,  z,  Uccpf,  Ux“Uccpf 
UXUCCPFl  Plots  streamwise  velocity  pertubations 

user  input 

•  PFxxx  data  file 

•  number  and  spacing  of  contour  levels 

output 

•  contour  plot  with  level  table 

VORTICITY  Computes  the  three  mean  vorticity  components 
user  input 

•  Vxxx  data  file 
output 

•  Wxxx  data  file  with  y,  z,  o>e,  o>r,  uz 
VORTICITY 1  Plots  any  desired  vorticity  contour 

user  input 

•  Wxxx  data  file 

•  number  and  spacing  of  contour  levels 

•  minimum  threshold  intensity  to  plot 

•  type  of  vorticity  to  plot  (x,y  or  z) 

output 

•  contour  plot  with  level  table 
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3DPLOTlxxx  Plots  a  3-D  array  of  up  to  seven  of  the  following  type  of 
plots.  Urn  plots  streamwise  velocity.  Ucc  plots  velocity  pertubations. 
Pt  plots  Pamb  minus  Ptotal.  Vx  plots  Mg.  Vy  plots  o>r .  Vz  plots  oz.  All 
user  inputs  are  identical  to  those  listed  in  the  previously  mentioned 
plotting  programs  with  the  additional  information  needed,  i.e.,  the  number 
of  plots  and  the  identification  of  all  data  files  to  plot.  Output  is  a  3- 
D  plot  of  all  graphs. 

C.  12  ELEMENT  HOT-FILM  PROBE  PROGRAMS 

12FILM1  Acquires  raw  and  conditioned  bridge  output  voltage  for  six 
active  hot-film  probes. 
user  input 

•  active  probe  numbers 

•  De 

output 

•  overall  sampling  frequency 

•  individual  channel  sampling  frequency 

•  FDATAxxx  data  file  with  active  probes,  raw  anc  conditioned  voltages 
C0V2  Analyzes  data  acquired  by  12FILM1 

user  incut 

•  FDATAxxx  data  file 
output 

•  mean  raw  voltages 

•  conditioned  rms  voltages 

•  normalized  rms  voltages 

•  covariances 

•  correlations 

•  correlation  versus  z/d  plots  for  increasing  and  decreasing  z 

D.  HOT-WIRE  PROBE  PROGRAMS 

PROBE003  Calibration  data  for  the  hot-wire  probe  containing:  Tcal, 
A,  ohr,  n,  B,  EOC,  Uljm  and  coefficients  AO,  Al,  A2,  and  A3 

CNT_BUF_AD  Data  file  used  to  configure  the  HP  6944A  multiprogrammer 
supplied  with  the  instrument 
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HOTWIRESPM  Computes  instantaneous  streamwise  velocities  and  velocity 
fluctuations.  Stores  5000  instantaneous  velocities  in  a  data  file. 
user  input 

•  sampling  frequency 

•  barometric  pressure 

•  ambient  temperature 

•  bulk  velocity 

•  bridge  output  zero  offset 

•  signal  conditioner  gain 

•  PRQBB003  data  file 

•  analog  mean  voltage 

•  analog  rms  voltage 

output 

•  WDATxxxxx  data  file  containing  5000  Ufl 

•  channel  Re 

•  De 

•  digital  average  Un 

•  digital  u'r  6 

HOTWIREVEL  Computes  normalized  velocities  both  digital  and  analog. 
Plots  U0  versus  time. 
user  input 

•  WDATxxxxx  data  file 

•  U  bulk  for  the  given  De 

output 

•  analog  U 

•  analog  u'2 

•  analog  and  digital  u' 

•  analog  and  digital  u'/U_ 

•  analog  and  digital  u'2/U2 

•  plot  of  Ug 
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APPENDIX  C  -  DYNAMIC  SIGNAL  ANALYZER  OPERATION 


The  following  sections  list  sequential  operations  to  acquire  frequency 
response  and  power  spectra  plots  using  the  Hewlett  Packard  3562A  Dynamic 
Signal  Analyzer.  Keys  with  []  codes  are  permanent  hard  keys.  Keys  with 
{>  codes  are  user  input  soft  keys  with  CRT  menus.  Keys  with  ()  codes  are 
permanent  hard  keys  shifted  to  alpha-numeric  code. 


A.  FREQUENCY  RESPONSE 

A  power  spectrum  plot  from  frequency  response  is  obtained  as  follows: 

•  Power  up  the  dynamic  Bignal  analyzer 

•  Connect  the  input  BNC  connector 

•  [MEAS  MODE] 

•  {LINEAR  RES} 

•  ( SELECT  MEAS ] 

•  {FREQ  RESP} 

•  [MEAS  DISPLAY] 

•  {POWER  SPEC1} 

•  [FREQ] 

•  {FREQ  SPAN} 

•  [0][.][1][,][2][5] 

•  {Hz} 

This  sets  up  a  frequency  response  from  0.1  to  25  Hz. 

•  [ SOURCE ] 

•  {SOURCE  OFF} 

•  [WINDOW] 

•  {HANN} 

•  [AVG] 

•  {STABLE  (MEAN)} 

•  {NUMBER  AVGS} 

•  [5] 

This  sets  up  5  stable  averages. 

•  [A] 

•  [SINGLE] 

•  [COORD] 

•  {MAG  (LOG)} 

•  {NEXT} 

•  {LOG  X} 

•  {RETURN} 

•  [UNITS] 

•  {P  SPEC  UNITS} 

•  {V2/Hz  (PSD)} 

•  {RETURN} 

•  {Hz  (SEC)} 

•  {TRACE  TITLE} 
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(D) (E) (A) (N) (BL) (=) (BL) [1] [0] [0] 
{ENTER} 


This  makes  the  title  DEAN  =  100. 


•  [SCALE] 

•  {X  FIXD  SCALE} 

•  [0][.][1][,][2][5] 

•  {Y  FIXD  SCALE] 

•  [0][.][0][0][0][0][0][0][0][0][0][0][1][,][1][0][0] 

•  [STATE /TRACE] 

This  checks  the  setup. 

•  [STATE/TRACE] 

•  [START] 


Wait  until  the  CRT  indicates  the  measurement 


is  complete. 


Ensure  the  plotter  is  ready. 


•  [PLOT] 


B.  POWER  SPECTRUM 

A  power  spectrum  plot  is  obtained  as  follows: 

•  Power  up  the  dynamic  signal  analyzer 

•  Connect  the  input  BNC  connector 

•  [MEAS  MODE] 

•  {LINEAR  RES} 

•  [SELECT  MEAS] 

•  {POWER  SPEC} 

•  {CHI  ACTIVE} 

•  [MEAS  DISPLAY] 

•  {POWER  SPEC1 } 

•  [FREQ] 

•  {FREQ  SPAN} 

•  [0] [ . ] [0] [0] [0] [1] [ , ] [1] [0] 

•  {kHz} 

This  sets  up  a  frequency  response  from  0.1  to  10000  Hz. 

•  [ SOURCE ] 

•  {SOURCE  OFF} 

•  [WINDOW] 

•  {HANN} 

•  [AVG] 

•  {STABLE  (MEAN)} 

•  {NUMBER  AVGS} 

•  (5] 

This  sets  up  5  stable  averages. 

•  [A] 

•  [SINGLE] 

•  [COORD] 

•  {MAG  (LOG)} 

•  {NEXT} 

•  {LOG  X} 
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•  {RETURN} 

•  [UNITS] 

•  {P  SPEC  UNITS} 

•  {V2/Hz  (PSD)} 

•  {RETURN} 

•  {Hz  (SEC)} 

•  {TRACE  TITLE} 

•  (D) (E) (A) (N) (BL) {=) (BL) [1] [0] [0] 

•  {ENTER} 

This  makes  the  title  DEAN  =  100. 

•  [SCALE] 

•  {X  FIXD  SCALE} 

•  [0]M[1]M[1J[0][0][0][0] 

•  {Y  FIXD  SCALE] 

•  [0][.][0][0][0][0][0][0][0][0][0][0][1][,][1][0][0] 

•  [STATE/TRACE] 

This  checks  the  setup. 

•  [STATE /TRACE] 

•  [START] 

Wait  until  the  CRT  indicates  the  measurement  is  complete. 
Ensure  the  plotter  is  ready. 

•  [PLOT] 
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APPENDIX  C  -  TABLES 


Table  1.  12  ELEMENT  HOT-FILM  PROBE  LOCATIONS 


PROBE  # 

z/d 

REGION  CHARACTERISTIC 

1 

4.48 

POSITIVE  VORTICITY 

2 

4.73 

DOWNWASH 

3 

4.98 

NEGATIVE  VORTICITY 

4 

5.23 

UPWASH 

5 

5.48 

POSITIVE  VORTICITY 

6 

5.73 

DOWNWASH 

7 

5.98 

NEGATIVE  VORTICITY 

8 

6.23 

UPWASH 

9 

6.48 

POSITIVE  VORTICITY 

10 

6.73 

POSITIVE  VORTICITY 

11 

6.98 

DOWNWASH 

12 

7.23 

NEGATIVE  VORTICITY 
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Table  2.  VORTEX  PAIR  APPEARANCE  AND  DISAPPEARANCE  OCCURRENCES 
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Table  3.  VORTEX  PAIR  APPEARANCE  AND  DISAPPEARANCE  FREQUENCIES 


APPEARANCES 

DISAPPEARANCES 

Dean 

e 

Grow 

Split 

Total 

Collapse 

Merge 

Engulf 

Cancel 

Total 

60 

85 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

O.CO 

0.00 

95 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

105 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

115 

0.74 

0.00 

0.74 

0.74 

0.00 

0.00 

0.00 

0.74 

125 

0.08 

0.00 

0.08 

0.08 

0.00 

0.00 

0.00 

0.08 

135 

0.98 

0.00 

0.98 

0.90 

0.02 

0.00 

0.04 

0.96 

75 

85 

0.47 

0.00 

0.47 

0.47 

0.00 

0.00 

0.00 

0.47 

95 

1.12 

0.00 

1.12 

1.13 

0.00 

0.00 

0.03 

1.16 

105 

1.49 

0.02 

1.51 

1.36 

0.02 

0.00 

0.00 

1.38 

115 

0.58 

0.14 

0.72 

0.33 

0.42 

0.00 

0.05 

0.80 

125 

1.70 

0.02 

1.72 

0.90 

0.30 

0.12 

0.40 

1.72 

135 

1.10 

0.02 

1.12 

0.57 

0.35 

0.09 

0.21 

1.22 

100 

85 

2.96 

0.05 

3.01 

2.30 

0.44 

0.00 

0.02 

2.76 

95 

2.31 

iltl 

2.48 

1.97 

0.56 

0.05 

0.15 

2.73 

105 

2.02 

0.10 

2.12 

0.73 

0.51 

0.25 

0.66 

2.15 

115 

2.03 

0.15 

2.18 

0.60 

0.60 

0.43 

0.42 

2.05 

125 

2.02 

0.00 

2.02 

0.42 

0.56 

0.51 

0.42 

1.91 

12o 

85 

1.07 

0.10 

1.17 

0,40 

0.33 

0.22 

0.09 

1.04 

95 

3.66 

0.16 

3.82 

1.45 

1.10 

0.52 

0.41 

3.48 

105 

1.38 

0.08 

1.46 

0.33 

0.33 

0.38 

0.18 

1.22 

150 

85 

2.48 

0.15 

2.63 

1.07 

0.76 

0.71 

0.17 

2.71 

95 

2.39 

0.21 

2.60 

1.20 

0.53 

0.89 

0.33 

2.95 
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Table  4.  VORTEX  PAIR  APPEARANCE 
FREQUENCIES 


AND  DISAPPEARANCE  NON-DIMENSIONAL 


APPEARANCES 


DISAPPEARANCES 


e 

Grow 

Spl  i  t 

Total 

Collapse 

Merge 

Engulf  Cancel 

Total 

85 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00  0.00 

0.00 

95 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00  0.00 

0.00 

105 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00  0.00 

0.00 

115 

.057 

0.00 

.057 

.057 

0.00 

0.00  0.00 

.057 

125 

.006 

0.00 

.006 

.006 

0.00 

0.00  0.00 

.006 

135 

.076 

0.00 

.076 

.070 

.002 

0.00  .003 

.075 

85 

.029 

0.00 

.029 

.029 

0.00 

0.00  0.00 

.029 

95 

.069 

0.00 

.069 

.070 

0.00 

0.00  .002 

.072 

105 

.093 

.001 

.094 

.084 

.001 

0.00  0.00 

.085 

115 

.036 

.009 

.045 

.020 

.026 

0.00  .003 

.049 

125 

.105 

.001 

.106 

.056 

.019 

.007  .025 

.107 

135 

.068 

.001 

.069 

.035 

.021 

.005  .013 

.074 

85 

.138 

.003 

.141 

.107 

.020 

0.00  .001 

.128 

95 

.107 

.008 

.115 

.091 

.026 

.002  .007 

.126 

105 

.094 

.005 

.099 

.034 

.024 

,012  .031 

.101 

115 

.095 

.007 

.102 

.028 

.028 

.020  .019 

.095 

125 

.094 

0.00 

.094 

.019 

.026 

.024  .019 

.088 

85 

.004 

.044 

.015 

.012 

.008  .003 

.038 

95 

.136 

.006 

.142 

.054 

.041 

.019  .015 

.129 

105 

.052 

.003 

.055 

.012 

.012 

.014  .007 

.045 

85 

.077 

.005 

.082 

.033 

.024 

.022  .005 

.084 

95 

.074 

.007 

.081 

.037 

.016 

.027  .010 

.090 

Table  5.  HOT-WIRE  PROBE  SENSING  LOCATIONS  RELATIVE  TO 


VORTEX  PAIR  STRUCTURE 


PROBE  POSITION 

z/d 

y/d 

REGION  CHARACTERISTIC 

1 

6.30 

0.10 

UPWASH  BASE 

2 

6.30 

0.50 

UPWASH  CENTER 

3 

6.95 

0.10 

DOWNWASH  BASE 

4 

6.60 

0.10 

BENEATH  VORTEX 

5 

6.95 

0.50 

DOWNWASH  CENTER 
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Table  6.  HOT-WIRE  SPECTRA  SIGNAL  CONDITIONER  SETTINGS 


Type 

meas. 

Dean 

SIGNAL  CONDITIONING  UNITS 

1.  Analog  E 

2.  Analog  e'  | 

Filter 

Gain 

Filter 

Gain 

High 

Low 

High 

Low 

Data 

acg. 

<275 

DC 

10  kHz 

2.0 

0.1  Hz 

1.0  kHz 

500* 

>275 

DC 

3  kHz 

2.0 

LF  plot 

All 

0.1  Hz 

0.1  kHz 

100 

HF  plot 

0.1  Hz 

3  kHz 

100 

Or  as  high  as  possible  to  not  saturate  the  rms  unit 
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APPENDIX  E  -  FIGURES 


Figure  1.  Transparent  Curved  Channel 
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Figure  2.  Blower  Assembly 
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SCREENS 


PIENCM  tZ 


Figure  5.  Five-hole  Pressure  Probe 


Figure  7.  Five-hole  Pressure  Probe  Dimensions  and  Connections 
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Figure  8.  Five-hole  Survey  Data  Acquisition  and  Traversing  Control 
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Figure  11.  Constant  Temperature  Anemometry  Bridge  and  Low  Speed  Data 
Acquisition  Systems 
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Figure  23.  Velocity  Perturbation  Contours  De=50.4  to  De=200.3 
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Figure  25.  Velocity  Perturbation  Contours  De=275.2  to  Pe=424.8 
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Figure  26.  Streamwise  Vorticity  Contours  De=50.4  to  De=200.3 
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Figure  29.  Radial  Vorticity  Contou.s  De=50.4  to  De=200.3 
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Figure  3C.  Radial  Vorticity  Contours  De=175.2  to  De=325.5 
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Figure  31.  Radial  Vorticity  Contours  De=275.2  to  De=424.8 
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Figure  36.  Correlation  vs.  z/d  Probes  2,4, 6, 8, 10, 12, wrt  2,  De«50.25 
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Figure  37.  Correlation  vs.  s/d  Probes  2,4, 6, 8, 10,12,wrt  2,  De«100 
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Figure  39.  Correlation  vs.  z/d  Probes  2,4,6,8,10,12,wrt  2,  De=200.7 
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Figure  44.  3-D  Correlations  vs.  z/d 
De=50.25  to  De=399.8 
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Figure  53.  3-D  Correlations  vs.  z/d  Probes  2,4,6,8,10,12,wrt  12,  from 
De=50.25  to  De=399.8 
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Figure  54.  Average  Vrms/Vmean  vs.  De  Probes  2,4,6,8,10,12 
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Figure  55.  Vrms/Vmean  /  VrmB399.8/Vmean399.8  vb.  r/d  from  De*50.25  to 
De=399.8 
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Figure  74.  3-D  Correlations  vs.  z/d  Probes  1,2, 3, 4, 5, 6,  wrt  6,  from 
De=49.2  to  De=400.3 
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Figure  76.  Vrm8/Vmean  /  VrtnB400.3/Vmean400.3  vs.  z/d  from  De=49.2  to 
De=400. 3 
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Figure  77.  Length  Scales  vs.  De  Probes  1,2, 3, 4, 5, 6 
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Figure  81.  Analog  Longitudinal  Turbulence  Intensity  vs.  De, 
Positions  1,2, 3, 4, 5 
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Figure  83.  Analog  Normalized  Longitudinal  Turbulence  Intensity 
vs.  De,  Positions  1,2, 3, 4, 5 
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Figure  84.  Hot-wire  Probe  Low  Frequency  Power  Spectra  Settings 
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Figure  152.  Hot-wire  Probe  High  Frequency  Power  Spectra  Settings 
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Figure  215.  Vortex  Pair  Split  Visualization,  De=100,  0=115° 
Interval  1/30  s 
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Figure  219. 
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Vortex  Pair  Appearance  and  Disappearance  Frequencies  vs.  De, 
6=85° 


147 


ADNanbani  ivNoisNawa  non 


Figure  220.  Vortex  Pair  Appearance  and  Disappearance  Non-dimensional 
Frequencies  (irfd/U)  vs.  De,  0=85° 
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Figure  221.  Vortex  Pair  Appearance  and  Disappearance  Frequencies  vs.  De, 
|  6=95° 
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Figure  224.  Vortex  Pair  Appearance  and  Disappearance  Non-dimensional 
Frequencies  (nfd/U)  vs.  De,  0=105° 
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Figure  226.  Vortex  Pair  Appearance  and  Disappearance  Frequencies  vs.  0, 
De=100 
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Figure  227.  Vortex  Pair  Appearance  and  Disappearance  Frequencies  vs. 
0=  85°  to  135°  Composite  w/o  Curve  Identification 
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Figure  228. 
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Figure  '>'>9.  Vortex  Pair  Appearance  and  Disappearance  Frequencies  vs.  De, 
6=  115°  to  135° 
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Figure  230.  Vortex  Pair  Appearance  and  Disappearance  Non-dimensional 
Frequencies  (nfd/U)  vs.  De,  0=85°  to  135°  Composite  w/o 
Curve  Identification 
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Figure  231.  Vortex  Pair  Appearance  and  Disappearance  Non-dimensional 
Frequencies  (trfd/0)  vs.  De,  0=85°  to  105° 
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Figure  232.  Vortex  Pair  Appearance  and  Disappearance  Non-dimensional 
Frequencies  (wfd/U)  vs.  De,  0=115°  to  135° 
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Figure  92.  Low  Frequency  Power  Spectrum,  De=50,  Position  1 
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Figure  93.  Low  Frequency  Power  Spectrum,  De=100,  Position  1 
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Figure  100.  Low  Frequency  Power  Spectrum,  De=300,  Position  1 
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Figure  103.  Low  Frequency  Power  Spectrum,  De=450,  Position  1 
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Figure  105.  Low  Frequency  Power  Spectrum,  De=lOO,  Position  2 
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Figure  110.  Low  Frequency  Power  Spectrum,  De=200,  Position  2 
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Figure  111.  Low  Frequency  Power  Spectrum,  De=250,  Position  2 
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igure  113.  Low  Frequency  Power  Spectrum,  De=350,  Position  2 
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Figure  115.  Low  Frequency  Power  Spectrum,  De=450,  Position  2 
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Figure  123.  Low  Frequency  Power  Spectrum,  De=250,  Position  3 
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Figure  126.  Low  Frequency  Power  Spectrum,  De=400,  Position  3 
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Figure  128.  Low  Frequency  Power  Spectrum,  De=50,  Position  4 
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Figure  132.  Low  Frequency  Power  Spectrum,  De=170,  Position  4 
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Figure  134.  Low  Frequency  Power  Sprctrum,  De=200,  Position  4 
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Figure  136.  Low  Frequency  Power  Spectrum,  De=300,  Position  4 
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Figure  141.  Low  Frequency  Power  Spectrum,  De=100,  Position  5 


210 


□- 

0. 


G)  01 
0  0 
J  2 


D  \ 
©CM 
0.  > 


□ 


D 

X 

lx. 


Figure  142.  Low  Frequency  Power  Spectrum,  De-150,  Position  5 
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Low  Frequency  Power  Spectrum,  De=160, 


Position  5 
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Figure  145.  Low  Frequency  Power  Spectrum,  De=180,  Position  5 
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Figure  146.  Low  Frequency  Power  Spectrum,  De=200,  Position  5 
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Figure  148.  Low  Frequency  Power  Spectrum,  De=300,  Position  5 
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Figure  157.  High  Frequency  Power  Spectrum,  De=150,  Position  1 
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Figure  159.  High  Frequency  Power  Spectrum,  De=170,  Position  1 
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Figure  163.  High  Frequency  Power  Spectrum,  De=300,  Position  1 
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Figure  177.  High  Frequency  Power  Spectrum,  De=400,  Position  2 
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Figure  183 


.  High  Frequency  Power  Spectrum,  De=170, 


Position  3 
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Figure  184.  High  Frequency  Power  Spectrum,  De=180,  Position  3 
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Figure  185.  High  Frequency  Power  Spectrum,  De=200,  Position  3 


251 


FxdXY  □  Log  Hz:  DEAN 


Figure  197.  High  Frequency  Power  Spectrum,  De=200,  Position 
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Figure  199.  High  Frequency  Power  Spectrum,  De=300,  Position  4 
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Figure  200.  High  Frequency  Power  Spectrum,  De=350,  Position  4 
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Figure  203.  High  Frequmcy  Power  Spectrum,  De=50,  Position  5 
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Figure  205.  High  Frequency  Power  Spectrum,  De=150,  Position  5 
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Figure  207.  High  Frequency  Power  Spectrum,  De=170,  Position  5 
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Figure  208.  High  Frequency  Power  Spectrum,  De=180,  Position  5 
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